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PREFACE 


Le  calcul  des  dcoulements  tridimensionnels  instationnaires  transsoniques  sur  des  ailes 
supercritiques  est  un  travail  difficile  qui  requiert  pour  1 'instant  des  temps  d’ordinateur 
incompatibles  avec  une  exploitation  industrielie.  Messieurs  Couston,  AngtSlini  et  Meurzec 
ont  pr<5sent6  a  la  Sous-Commission  d’Adroelasticiti  du  Panel  une  approche,  fondle  sur  la 
m£thode  des  tranches,  dans  laquelle  I’aile  est  template  par  une  succession  de  profils  dont  les 
interactions  sont  prises  en  compte  d’une  maniere  approchfe;  les  temps  de  calculs  deviennent 
alors  raisonables, 

Les  auteurs  ont  prdsente1  (’application  de  la  methode  au  calcul  des  champs  de  pression 
instationnaire  induits  par  1’oscillation  de  I'aileron  de  l'aile  ZKP,  dans  les  conditions  de  1'essai 
effectue  sur  cette  aile  supercritique  dans  la  soufflerie  SI  Ma  de  Modane;  la  eomparaison 
theorie -experience  n’est  pas  parfaite,  du  fait  que  les  cffets  visqueux  n’ont  pas  6t£  pris  en 
compte;  elle  est  cependant  encourageante  dans  le  sens  ou  les  principaux  phenomenes  sont 
correctement  pr&lits. 

La  Sous-Commission  d’A£ro£ last icit<5  est  heureuse  de  publier  ce  rapport,  qui  est 
I’aboutissement  d’un  programme  cooperatif  patronniJ  par  le  Ministere  Allemand  de  la 
Recherche  et  Technologic  et  1’ONERA,  avec  la  participation  du  DFVLR,  de  MBB,  de  VFW- 
Fokker  et  de  1’ONERA. 


G.COUPRY 

Chairman,  Sub-Committee 
on  Aeroelasticity 
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COMPARISON  OES  CHAMPS  DE  PRESSION  INSTATIONNAIRES 
CALCULES  ET  MESURES  SUR  LE  MODELE  ZKP 

par  Michel  COUSTON,  JeanJacques  ANGELINI  et  Jean  Louis  MEURZEC 

Office  National  d' Etudes  et  de  Recherches  Airospatiaies  (ONER A) 
92320  ChStillon  (France) 


RESUME 


La  prediction  dee  champs  de  pression  instationnaire  sur  une  aile  supercritique  est  envisagde 
en  utilisant  une  mdthode  bidimensionnelle non  lindaire  corrigde  des  interactions  tridimensionnelles 
entre  tranches  A  partir  d'une  approximation  asymptot ique .  Cette  approximation  est  applicable  aux  ailes 
de  grand  allongement  et  de  fliche  moddrde. 

Une  partie  de  cet  article  porte  sur  une  application  de  la  mdthod e  £  un  cas  reel,  ce  cas 
correspond  aux  essais  dans  la  soufflerie  transsonique  Sl-Modane  d'un  module  rigide  dquipe  d'une 
gouverne  oscillante.  Ces  essais  font  partie  du  programme  ZKP  qui  a  ete  patronnd  a  la  fois  par  le 
Ministbre  de  la  Recherche  et  de  la  Technologie  de  la  Rdpublique  Federate  Allemande  et  l'ONERA  avec 
la  participation  du  DFVLR,  de  MBB,  de  VFV-FOKKER  et  de  l'ONERA. 

SUMMARY 


The  unsteady  pressure  fields  prediction  for  a  supercritical  wing  is  considered  from 
a  bidimensional  nonlinear  method  corrected  for  tridimensional  interactions  between  strips  using 
an  aaymptotic  approximation.  This  approximation  is  derived  for  wing  with  large  aspect  ratio  and 
low  sweep  angle. 

Part  of  the  paper  will  deal  with  an  application  of  this  method  to  a  real  case,  that  is 
the  unsteady  experiment  carried  out  with  oscillation  of  the  control  surface  on  a  rigid  model  on 
the  Sl-Modane  transonic  wind-tunnel.  These  experiments  have  been  part  of  the  ZKP  programme , 
sponsored  by  the  German  Federal  Ministry  of  Research  and  Technology  and  ONERA,  with  participation 
of  DFVLR,  MBB,  VFW-FOKKER  and  ONERA. 
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I  -  INTRODUCTION 


La  prediction  des  dcoulements  instat ionnaires  sur  les  ailes  aupercrit iques  dquipant  les 
aviona  de  transport  actuels  eat  un  probldme  part icul idrement  ddlicat.  La  presence  d'ondes  de  choc 
remet  en  cause  lea  hypotheses  traditionnellea  de  lindaritd  des  phPnotnines.  Pourtant,  que  l'on 
envisage  une  etude  de  stabilite  pour  prdvenir  un  eventuel  flottement  ou  que  l'on  s'interesse  aux 
techniques  de  Contrble  Actif  Generalise  (C.A.G  ou  C.C.V.  en  anglais),  il  faut  connaftre  les  charges 
aerodynamiques  instat ionnaires  agissant  sur  la  structure. 

Dans  cette  optique,  l'ONERA  a  entrepris  des  etudes  experimentales  sur  les  effets  instation- 
naires  lies  A  1 'oscillation  des  gouvernea.  Ces  etudes  d'abord  bidimens ionnel les  [l,2]  ont  ete 
etendues  A  une  aile,  dans  les  deux  caa  il  s'agisaait  d'un  profil  supercritique  developpe  par 
1 'Aerospatiale,  (type  RAI6SC1 ) .  Puis  en  collaboration  avec  VFW-FOKKER,  MBB ,  le  DFVLR  et  le 
MinistAre  Allemand  de  la  Recherche  et  de  la  Technolologie  les  champs  des  pressions  indationnaires 
ont  egalement  ete  mesurds  sur  le  models  ZKP  qui  est  l'une  des  conf igurations  seroeiastiques 
tridimensionnel  les  retenues  par  l'AGARD[3]  . 

Paral lAlement  A  ces  etudes  experimentales,  des  progrAs  interessants  ont  ete  obtenus  sur  le 
plan  thdorique  en  bidimenaionnel  [  4 , 5 ,6 J  .  Une  approche  tr idimensionnel le  non  1 indai re  est  en  cours 
de  ddveloppement  A  la  Direction  des  Structures  de  l'ONERA  mais  elle  ndcessiters  des  temps  de 
calcul  beaucoup  plus  important  qu'en  bidimenaionnel.  On  propose  done  ici  une  approche  provisotre 
basde  sur  l'idde  de  1 ' approximat ion  asymptotique  pour  les  ailes  de  grand  allongement  et  faible 
flAche  qui  permet  d'introduire  une  correction  tridimensionnel le  dans  le  calcul  bidimenaionnel 
non  lindaire. 

On  presente  ensuvte  une  compataison  entte  des  Tdsultats  experimentaux  et  des  calculs  non 
1 indaires  bidimens ionnels  et  tridimensionnels  approchds.  Cette  comparsison  s'appuie  principalement 
sur  l'aile  ZKP  esaayde  dans  la  Soufflerie  Sl-Modane  de  l'ONERA. 


II  -  METHODE  DE  CALCUL  DES  EfFETS  INSTAT IONNAIRES 
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II. 1  -  Mdthode  bidimenaionnelle  de  bate 

La  resolution  numdrique  dea  effets  inatationnairea  aur  une  aile  supercritique  ndceaaite 
un  calcul  tridimensionnel  non  lindaire.  Dea  efforta  aont  en  coura  4  la  Diviaion  d'Adrodlaaticitd 
de  l'ONERA  pour  ddvelopper  un  Cel  code  de  calcul  &  partir  dea  travaux  prdaentda  en  rdfdrence  [5,6] 
mais  ce  code  de  calcul  n'eac  paa  diaponible  dana  l'imraddiat  et  11  ndceaaitera  dea  temps  de  calcul 
assez  importanta .  On  ae  propoae  done  pour  1'inatanC  d'utiliaer  la  mdthode  bidimenaionnelle  [  6]  en 
lui  apporCant  dea  correctiona  tridimenaionnellea . 

La  mdthode  bidimenaionnelle  ne  sera  paa  ddcrice  en  ddtail  daracet  article  ayant  ddji  fait 
l'objet  de  publicationa  antdrieurea  [5,6]  .  Noua  rappelerona  qu'elle  permet  de  rdaoudre  l'dquation 
dea  petitea  perturbationa  tranaaoniquea  et  inatationnairea  (1)  : 


a* mi  tl  tZxmi  .  A|f  1  £»  + 

m  di'*  “3757  L  3«->37  3^ 


X  r[  (t  ♦  1  )  m  i  ♦  3(4 -••£)]  "L 


La  resolution  de  liquation  (la)  est  obtenue  par  une  extension  de  la  mdthode  implicite 
de  directions  alterndes  proposde  par  BALLHAUS  .  Le  schema  numdrique  se  decompose  en  trois  pas 
et  se  prdsente  sous  la  forme  suivante  : 


-  pas  en  x 


-  pas  en  Y 
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•  paa  de  calcul  de  la  variable  Z 


,"•*  -  y  " 

ZM - zd, 


+  is. 

J’k  +  at 


ou  n ,  j  ,  It  ,  aont  reapectivement  lea  indicea  du  temps,  de  I’abaciaae  et  de  1 'ordonnde  ;ou  le 
maillage  utiliad  eat  cartdaien  (lea  conditiona  aux  limites  du  profil  aont  dcritea  aur  un 
aegment  de  droite  reprdaentant  une  approximation  de  la  poafcion  moyenne  du  profil). 

Lea  aymbolea  fl«  et  dyy  reprdaentent  reapectivement  la  diffdrence  finie  ddcentrde 
en  arridre  auivant  «  et  la  ddrivde  aeconde  centrde  auivant  i 

D,„  eat  un  opdrateur  mixte  aimilaire  4  celui  de  MURMAN  [8]  qui  assure  la  diserdti- 
aation  du  terme  non  lindaire  en  prdaervant  4  la  foia  la  atabilitd  du  schdma  (par  un  ddeentrement 
dana  lea  zonea  auperaoniquea ) , la  conaervativitd  (par  un  traitement  spdeial  des  points  aoniques  et  des 
points  "choc")  et  la  lindarisation  du  probldme  (par  une  diaerdt isation  temporelle  appropride  du 
terme  non  lindaire  entre  lea  deux  tempa  n  et  (  n.  ♦  l)At  ). 

Notons  que  le  schdma  numdrique  (2a  et  2b)  est  intdressant  car  il  ne  gdn^re  que  dea 
matrices  tri  ou  quadridiagonalea  qui  aont  rdsolues  efficacement  par  double  triangulariaation. 

Le  troisidme  paa  (2c)  calcule  la  variable  Z  point  par  point  et  ne  ndeessite  pas  de  rdsolution 
matricielle .  Ce  troisidme  paa  a  dtd  introduit  pour  traiter  la  ddrivde  aeconde  temporelle  de 

manidre  centrde  ce  qui  n'est  paa  possible  dans  un  schdma  4  deux  pas  comne  celui  utiliad  en  [$]  . 

11.2  -  Champ  atationnaire 

Avant  de  faire  un  calcul  inatationnaire  il  faut  ndeessairement  avoir  un  champ  de  pression 
atationnaire  raiaonnable.  En  effet,  puiaque  1'on  conaiddre  ici  le  cas  d'ailes  opdrant  en  rdgime 
tranaaonique  ,  il  aerait  illuaoire  de  faire  un  calcul  inatationnaire  ai  la  position  du  ou  dea 
ondea  de  choc  n'eat  paa  correcte.  On  doit  done  prddire  de  manidre  correcte  la  position  dea 
zonea  auperaoniquea.  Pour  ce  faire,  il  faut  connaltre  pour  cheque  section  d'aile  ou  le  calcul  eat 

demandd  lea  paramdtres  4  fournir  au  calcul  bidimenaionnel .  Ces  paramdtrea  atat ionnairea  aont 

le  nombre  de  Mach,  1 'incidence  et  une  dventuelle  correction  pour  tenir  compte  de  la  ddportance 
induite  par  lea  effets  viaqueux  atat ionnairea . 

11.2.1  -  Correct ion_d^incidence_induite 

La  thdorie  de  I'aile  d'envergure  limitde  est  connue  depuia  lea  travaux  de  Prandtl,  elle 
nous  dit  que  ('incidence  effective  diffdre  de  l'incidence  adrodynamique  de  I'aile  4  cause  de 
la  vitease  induite  par  une  nappe  tourbi 1 lonnaire .  Cela  ae  traduit  dana  une  section  donnde  par  un 
angle  induit  qui  ddpend  de  Involution  (auivant  1’envergure)  de  la  portance.  On  doit  done  estimer 
la  valeur  de  cet  angle  induit  atationnaire. 
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Ce  calcul  va  Stre  fait  k  partir  de  la  mEthode  tridimensionnelle  des  doublets  ("doublet 
Lattice-method")  dEveloppEe  k  1 'ON ERA  k  partir  des  travaux  d 1  ALB ANO  et  RODEEN  [9]  .  Cette  mEthode 
ne  prend  pas  en  compte  les  effeta  d’Epaisseur,  done  si  1 ' on  veut  introduire  des  effets  de  portance 
il  faut  considErer  des  sections  d'aile  planes  mises  en  incidence. 

A  partir  du  vrillage  gEomEtrique  de  l'aile  ZKP  tel  qu'il  est  dEfini  en  [3]  et  repris 
en  figure  1,  nous  considErons  le  vrillage  de  plaques  planes  Equivalent  en  introduisant  1* angle  de 
portance  nulle  de  chaque  section  et  une  correction  empirique  qui  tient  compte  des  effets  visqueux. 
Ceci  nous  donne  la  loi  de  vrillage  prEsentEe  en  figure  1,  loi  qui  est  introduite  dans  le  programs 
linEaire  tridimens ionnel  des  doublets.  On  notera  que  l'aile  est  supposEe  infiniment  rigide  (d'autant 
qu’elle  est  haubanEe)  ce  qui  nous  autorise  k  nEgliger  le  dEvrillage  aErodynamique  (dans  le  cas  d'un 
avion, l'aile  n'est  pas  rigide  et  il  faudrait  tenir  compte  de  sa  dEformation  sous  la  charge). 

L'Evolution  du  coefficient  de  portance  ainsi  calculE  est  prEsentE  en  figure  2,  on  le 
compare  aux  mesures  obtenues  dans  Sl-Modane.  L' accord  est  relativement  sat isf aisant ,  nous  pouvons 
done  dEterminer  1' incidence  de  chaque  section  en  envergure  en  comparant  le  C|  bidimensionnel 
de  plaque  plane  au  C  (  tridimens ionnel  calculE  de  la  figure  2.  Nous  obtenons  la  loi  d* incidence 
sur  chaque  profil  que  nous  prEsentons  en  figure  3  (cette  incidence  est  dEfinie  par  rapport  k  la 
corde  de  chaque  section).  On  notera  que  la  section  du  bout  de  plume  ( T)  =  .885)  oil  sont  prElevEes 
les  pressions  instationnaires  travaille  k  une  incidence  faible  pour  la  configuration  choisie 
(  M  . ,  *  0.78,  tt  . ,  *  0®,  6  •,  *  0°),  nous  cons idererons  done  que  dans  cette  configuration 

la  section  travaille  sans  incidence. 

Il  est  clair  que  toute  variation  d' incidence  globale  de  l'aile  (ou  tout  braquage  de  la 
gouverne)  entraine  des  modifications  dans  l1 incidence  de  chaque  section  et  quels  courbe  prEsentEe 
en  figure  3  est  alors  modifiEe.  A  partir  du  programme  tridimens ionnel  linEaire  nous  dEterminons  ces 
variations,  ainsi  en  ne  considErant  que  la  section  c  .885  qui  est  instrumentE  avec  des  capteurs 
instationnaires  on  trouve  que  les  corrections  d' incidence  locale  prennent  la  forme  : 

-  pour  un  braquage  vers  le  haut  (nEgatif)  de  la  gouverne  d'un  degrE  1 ' incidence  locsle  pour 
cette  section  augmente  d' environ  0.395°. 

-  pour  une  mise  en  incidence  de  1* ensemble  de  l'aile  d'un  degrE  la  section  X\  *  .885  voit  son 
incidence  effective  croitre  seulement  de  0.725°. 

Ces  deux  derni&res  corrections  seront  appliquEes  par  la  suite  k  un  exemple  concret  compor- 
tant  k  la  fois  une  raise  en  incidence  de  l'aile  et  un  braquage  de  la  gouverne. 

I I . 2 . 2 . -  DEgor tance_ indui te_gar_ lesef fe tsvisqueux 

Il  s'agit  Ik  d’une  correction  empirique  qui  ne  se  justifie  qu'en  regard  d'une  comparaison 
avec  des  essais.  Nous  avons  observE  en  rEfErencetb]  que  pour  le  profil  bidimensionnel  supercritique 
RA16SC1  de  1 ' aErospatiale  la  portance  calculEe  Etait  systEmat iquement  surEvaluEe  en  stationnaire 
et  ce  k  cause  de  la  couche  limite  qui  se  dEveloppe  sur  le  profil.  Il  en  rEsulte  principalement  que 
les  ondes  de  choc  se  placent  tropen  arriEre  sur  le  profil;  pour  pallier  cet  inconvEnient  la  seule 
solution  qui  amEliore  la  comparaison  sans  nEcessiter  1 ' introduction  de  la  couche  limite  nous  semble 
etre  un  braquage  de  la  gouverne  vers  le  haut.  En  rEfErence[3»  nous  avons  montrE  qu'un  braquage  de 
-  2°  (vers  le  haut)  de  la  gouverne  donnait  de  bons  rEsultats  ;  ici  n'ayant  pas  d'essais  bidimension- 
nels  pour  les  profils  Equipant  l'aile  ZKP  nous  appliquerons  la  meroe  correction  puisque  les 
sections  sont  de  type  supercritique  tout  comme  le  RA16SC1  de  1 'AErospatiale  . 

I I. 2. 3  -  Correction  du  nombre  de  Mach 


Pour  faire  le  calcul  bidimensionnel  il  faut  se  donner  le  nombre  de  Mach  incident  pour 
la  section  considErEe.  En  rEfErence  flo]  on  trouve  certaines  reraarques  en  ce  qui  concerne  le 
nombre  de  Mach  k  considErer  pour  Evaluer  le  champ  de  pression  stationnaire  dans  la  section  droite 
d'une  aile  cylindrique  en  attaque  oblique  (d'un  angled  ).  Ce  nombre  de  Mach  est  thEoriquement  le 
Mach  de  l'aile  multipliE  par  co5»p  alors  que  des  travaux  expErimentaux  rEfErencefllJ  semblent 
indiquer  que  le  coefficient  serait  plutot  yco5<P  pour  une  aile  en  flEche.  En  suivant  les  remarques 
faites  sur  le  sujet  en  rEfErence  £lC^  et  en  appliquant  les  rEgles  de  ftandtl-Cl auert  pour  reus 
ramener  k  une  section  dans  le  lit  du  vent  nous  avons  obtenus  que  le  nombre  de  Mach  du  calcul 
bidimensionnel  Etait  donnE  par  la  formule  suivante  : 


Mw  -  .  I  >1  4  -wj.u  (3) 

V  coi1? 

L' Equation  (3)  ne  prEsente  aucun  problEme  pour  une  aile  cylindrique  ;  ici  pour  la 
maquette  ZKP  on  a  une  certaine  latitude  dans  le  choix  de  la  flEche  qui  n’est  pas  la  meme  au  bord 
d'attaque  et  au  bord  de  fuite.  Le  nombre  de  Mach  n'est  done  dEfini  qu'entre  deux  limites  et 
seule  I'expErience  peut  nous  permettre  de  trancher. 

I I. 2. 4  -  Aggl icat iondescorrect ionsstat ionnaires 

Nous  al Ions  considErer  deux  cas  expErimentaux  pour  valider  les  corrections  stat ionnaires 
telles  que  nous  venons  de  les  dEfinir,  ces  deux  essais  seront  repris  par  la  suite  en  instationnaire 
pour  une  excitation  narmonique  de  la  gouverne. 


X Is  correspondent  respectivement  h  : 


-  essai 

1  ; 

1  "  aile 

-  0.78 

;  a  ..  *o° 

aile 

;  6 

-  essai 

2  : 

:  "  aile 

-  0.78 

;  a  *  2° 

aile 

;  6 

Consid^rons  tout  d'abord  l'essai  1  .  Si  l'on  n' introduit  aacune  de6  corrections  stationnaires 

que  nous  venons  de  d^finir  on  obtient  le  r£sullat  de  calcul  prdsentd  en  figure  4  (pour  la  section 
H  *  .885  qui  est  la  section  instrument^  en  dynamique).  Nous  comparons  ce  r£sultat  aux  mesures 
dans  la  section  *  .85  (done  tr£s  voisine)  qui  elle,est  instruments  pour  les  pressions  statiques. 

On  observe  que  calcul  et  essai  sont  tr&s  diff£rents,  la  portance  calculde  est  ici  environ  quatre 
fois  trop  importante.  Si  nous  appliquons  la  correction  d£finie  pr£c£demment  sur  le  braquage  de  la 
gouverne  (2°  vers  le  haut)  et  que  nous  ndgligeons  la  correction  d' incidence  relativement  faible  pour 
cette  position  en  envergure  nous  obtenons  la  comparaison  prdsentde  en  figure  5.  L' accord  calcul- 
experience  est  meilleur,  on  observe  toutefois  une  onde  de  choc  de  part  et  d' autre  du  profil  dans 
le  calcul,  ondes  qui  ne  sont  pas  observdes  experiment a lenient .  Ces  ondes  de  choc  disparaissent  si 
l'on  introduit  la  correction  du  nornbre  de  Mach  en  tenant  compte  de  l'effet  de  flfcche,  au  nombre  de 
Mach  de  calcul  0.73  on  trouve  la  comparaison  presentee  en  figure  6.  L' accord  entre  les  mesures 

et  le  calcul  est  alors  suffisant  pour  envisager  1 'oscillation  harmonique  de  la  gouverne  autour 
de  ce  champ  moyen. 

Si  nous  envisageons  maintenant  l'essai  2  et  que  nous  appliquons  1 'ensemble  des  corrections 
stat ionnaires  sur  1' incidence,  le  braquage  de  la  gouverne  et  le  nombre  de  Mach  nous  obtenons  les 
donndes  suivantes  pour  le  calcul  bidimensionnel  : 

0.73  ;  =  3.82°  ;  &m  -  -  8° 


Ce  calcul  est  compard  aux  essais  en  figure  7.  On  voit  qu'ici  encore  le  champ  de  pression 
stationnaire  calculi  est  en  bon  accord  avec  le  champ  mesur£.  La  position  de  l'onde  de  choc  et  son 
intensity  sont  bien  rendus,  on  peut  done  dans  ce  cas  £galement  faire  osciller  la  gouverne  autour  de 
sa  position  moyenne  pour  g£n4rer  un  champ  de  pression  ins tationnaire ,  mais  auparavant  nous  devons 
introduire  des  corrections  instationnaires  tridimensionnelles . 

II. 3  -  Prise  en  compte  des  effets  tridimensionnels  lin^aires 

La  prise  en  compte  des  effets  3-D  dans  le  calcul  2-D  par  tranche  d'une  aile  de  grand 
allongement  s'avfcre  n£cessaire,  en  effet  on  sait  que  les  effets  instationnaires  sont  nettement 
accentuds  par  un  calcul  bidimensionnel.  On  peut  penser  que  la  difference  principale  entre  les  deux 
calculs  provient,  comme  le  proposait  Prandtl,  d'une  vitesse  ou  d'un  angle  induit  par  une  nappe 
tourbillonnaire ,  mais  ici  cet  angle  induit  serait  ins tationnaire .  Cet  effet,  purement  lindaire, 
existe  merae  si  1 ' aile  est  localement  supercritique  et  nous  allons  voir  qu'il  est  possible  d'estimer 
1 'angle  inuit  h  partir  des  m£thodes  lin£aires  existantes  et  de  corriger  de  manifcre  significative 
le  calcul  bidimensionnel  non  lin£aire. 

La  formulation  int£grale  du  problferae  lin£aire  par  la  m£thode  des  doublets  peut  s'£crire, 
en  utilisant  le  noyau  tridimens ionnel  D  ,  sous  la  forme  suivante  : 

,p3(x>-y<0)  =  +  = 

<4i 


Schema  1  :  Notations 


or  le  noyau  3-D  s'dcrit  : 


"30  dt 


(5) 


fw) 

est  ddfinie  &  l'aide  du  potentiel  d' acceleration  d'une  source  pulaanCe.  On  a  done  : 


=*“»*■  0  (  (  ACpU.v)  |L  ( 

>■’  “lb  J-!/t  3 3  J-co 


(6) 


soit  : 


►b 


3l  j  (  ACp  (u,,ir)  [  f3°(t/^-v,  j)  di.  d«dir 
J-.0  j  b  J_t/t  J-<*» 


(7) 


Le  probl&me  initial  est  la  resolution  de  l’dquation  lindaire  du  potentiel,  il  se  presente 
done  sous  la  forme  suivante  : 


(8) 


ou  sex  est  un  opdrateur  ne  dependant  que  de  la  variable  X  ,  en  ddrivant  suivant  ig 

(9) 

i  0 

1 

On  peut  done  dcrire  : 


n4*/l 

ACp(u.(v)  ?S0(t,^-V;2)  dtdu  cUr 


(10) 


Considdrons  tout  d’abord  le  premier  terme  du  second  membre  et  supposons  que  l’aile  est  de 
grand  allongement.  On  peut  alors  considdrer  que  le  ACp  varie  peu  en  envergure,  ce  qui  nous  permet 
d*dcrire  : 


,4  .X-JK.tCO 

=  '  &"»  3^x1  aCp(4i.,^)l  [  ^dtdudir 

>-?/j  J -00  J- to 


(11) 


ce  qui  ae  transforme  en  : 


♦?/,  x-M. 

=  3^x  ACp(at,^)  \  f*°(  l,l)  dtdu. 

3”*°  •'-*/»  J-oo 


(12) 


.4  tit  .Ot-M. 

«d(x,^)=  -  BLm  1  ACp(M,‘3)3fx  ?10(t,5)  dtdM. 

3~>0  J.t/j  J-co 

Ce  qui  si  lfon  consid&re  le  probldtne  bidimensionnel  lin£aire  : 

(3tx+|L)^°  =  ° 


nous  permet  d'dcrire  : 


et  done  : 


4  til  .X--U. 

L>»o  j^Cpf-u^)  |L  j  ■f?*°(  t, 


w*  (t 


*a)-  •v20(*--*>v  djA- 


(13) 


(14) 


(15) 


(16) 


Nous  reconnaissons  Id  la  formulation  intdgrale  du  probldme  lindaire  bidimensionnel , 
l'dquation  (10)  devient  alors  : 


I 


1-7 


+1/1 

w(x,^)=  lim  f  ACp(-u(-ij). 

/-hi. 


/-i/l 

i  (*iU 

-  —  I  \  ACp(j»,v)l  ?3°(  t  <m-v  alt  cl*  oLr 

3-°^llbJ-«4  .  J—  3  3 

Nous  pouvons  done  introduire  la  notion  d'un  angle  induit  qui  est  repr£sent£  par  le 

deuxi£me  terrae  du  second  membre  de  (17)  : 

+  S>  (  +  UL  /*-A 

=  —  I  (  1  f3°(t,4i.tf  o)  cLl  dm  djj  (18) 

*8  J-v,  J. e/i  0 

mais  si  l'aile  est  de  grand  allongement  on  peut  consid^rer  que  : 


|  ^-v  ;  °)  <°  j  (  bi  ‘#-v,  °)  °Lt 


car  -  fc/2,  <  Oc  -  M.  4  6/2. 


+  b  i*+  i/i.  ,  o 

=  —  I  &Cp  ( M,  v)  Am.  -f  (  fc  , 'VL-tf  o)  dt  djS  (20) 

^  J.b  J-e/i  J_.  d 


on  d£finit  : 


on  a  alors  : 


^30  =  1  ACp  (  ju.,  v)  dLu. 

J.  e/i 

(°  io 

K  M  =  i_  f*°(t,v)  dt 

3  V  J—  CO 

rS  : 

6i<^:|blir  ^30  K^-wjdu- 


On  voit  done  maintenant  que  Si  n'est  fonction  que  de  seul;Si  reprfsente  la 
correction  a  apporter  au  bidimensionnel  pour  prendre  en  compte  lea  effets  tridimensionnels ,  (17) 
devient  : 

,  +  tll 

'*'(*,'3)=  I  ^ )  cLu.  -  ( *3 )  (24> 

5~’°  -t'i 

Nous  pouvons  multiplier  (24)  par  une  fonction  4,(*)  dgfinir)  et  int£grer  le  long 
de  la  corde,  nous  obtenons  : 

/+e/i  /+e/t 


=  k  J _((1  I .^Cp(»,?iNlV*, 5) +**)<»’ 

,+tn. 

-  6i f^*)  J  ^  ♦  (x)  clx 


ou  en  changeant  X  en  -X  dans  le  second  membre  de  (25)  : 

C*tX  /♦***  C*tl  zo 

n"V-*,3)  *  (-x)d-dx  (26) 

,  +  iii 

-  J  f/  *(-*■)  <*•*• 

Si  nous  d£finissons  ^  (-3t)  comme  le  saut  de  pression  engendr£  par  un  pompage  nous 
obtenons  immldiatement  : 

.Air.  ,  +  l'i- 


f4  1  C 

W*(*)wU'^  °Lx  =  J.l<,aCP(4‘i*)  '  4*  Au- 


*< 


1-8 


y 


f- 


ou  est  la  pv.rtance  bidimensionnel le  attach^e  au  pompage.  Nous  obtenons  finalement  1 'Equation 

int£gro-di  f  ft§rentielle  (Prandtl  en  s tat ionnaire ,  Reissner  en  instationnaire)  : 


♦  bo(^)-  ^  —  fejo  K('M-'U)  ^  (28 

.{/j_  b  75 v 


qui  pennet  de  calculer  « 30  Puis  1‘ angle  induit.  Mais  disposant 
lineaire,  nous  pouvons  obtenir  la  portance  ^"50  directement 

(28)  la  valeur 

1  ft  ^30  -  j_ 

ft* 


tit 

t^4>(x)  w  (  a,  *j)  d.=c 


d’un  programme  tridimensionnel 
,  nous  tirons  done  de  l ’Equation 


(29) 


(30) 


(31) 


Done  en  ddfinissant  feio  comme  la  portance  bidimens ionnel le  de 
nous  obtenons  : 

»i(^=  jfc[  feao(<a)] 

Cette  derni£re  expression  amdne  quelques  remarques.  Tout  d'abord 
a  ,  e'est  done  un  pompage  dont  1* amplitude  c*  est  la  valeur  entre 

£galement  que  si  le  probl£me  tridimensionnel  est  quasiment  bidimensionnel 
(  r  ).  Si  est  1 'amplitude  du  pompage  on  a  : 


la  tranche  consid£r£e 


(32) 

est  proport ionnel 
crochets  ;  on  notera 
1' angle  induit  s'annule 


^•30  =  ^ 


20 


ftc 


(33) 


Cette  £galitd  indique  que  le  pompage  que  1 ' on  superpose  au  probl^me  bidimensionnel  sert 
h  r^tablir  la  valeur  correcte  de  la  portance.  On  peut  g£n£raliser  cette  id^e  en  prenant  pour 
valeur  au  lieu  d’une  fonction  constante  par  tranche  une  loi  lineaire  suivant  la  corde . 

Toute  loi  lineaire  suivant  la  corde  est  la  superpositon  d'un  pompage  et  d'un  tangage  ;  on 
d£terminera  done  les  deux  paramdtres  c*  et  ^  (qui  repr£sente  les  intensity  respectives 
du  pompage  et  du  tangage)  en  pr£servant  non  seulement  la  portance  mais  aussi  le  moment  au 
quart  avant  dans  la  section  soit  : 


ftjo=  ft  to  +  <*  ftet  (3  fcvj 

'mjo-  'W-to  +  +  P'm-b 


(34) 


et  Ttlb  sont  respect ivement  les  portance  et  moment  (quart  avant) 


li£s  au  tangage. 


La  resolution  du  syst£me  (34)  pour  determiner  les  valeurs  <*  et  {$  necesate  done  la 
connaissance  des  portances  et  moments  du  champ  tridimensionnel  lineaire,  ce  champ  est  calcuie 
h  partir  de  la  methode  tridimens ionnel le  instationnaire  dite  methode  des  doublets.  Connaissant 
ces  valeurs  on  injecte  dans  la  methode  non  lineaire  un  mode  d ’ osc i 1 lat ion  qui  est  la  superposition 
du  mouveraent  consider,  d'un  mode  de  tangage  au  quart  avant  et  d’un  mode  de  pompage  ;  ces  trois 
mouvements  etant  dephases  les  uns  par  rapport  aux  autres  puisque  o<  et  ^  sont  des  nombres 
complexes . 


Nous  allons  voir  dans  le  paragraphe  suivant  des  applications  de  cette  technique. 
Ill  -  COMPARAISON  DU  CALCUL  AUX  MESURES  INSTATIONNAI RES  SUR  L'AILE  ZKP 


Le9  me9ures  instationnaires  preievees  sur  la  maquette  ZKP  ont  ete  obtenues  pour  deux 
configurations.  La  premiere  correspond  h  l'aile  haubanee  h  son  extremite  avec  une  gouverne 
oscillant  de  mani&re  harmonique,  alors  que  pour  la  seconde  l'aile  est  libre  et  excitee  h  la  frequen¬ 
ce  d'un  mode  propre  par  1 'oscillation  de  la  gouverne.  Nous  limiterons  notre  eomparaison  h  la 
premiere  con f iguration  car  nous  ne  connaissons  pas  les  valeurs  des  filches  engendrees  par  le 
mode  propre.  Ce  faisant,  la  seule  section  de  mesures  instationnaires  oil  1 ' on  observe  des  evolutions 
de  pression  signi f icat ives  est  la  section  T)  *  .885.  Les  autres  sections  de  mesures  T)  *  .41 
et  T|  «  .64  sont  trop  eioignees  de  2a  gouverne  qui  avec  une  profondeur  de  22.5  7,  s'etend  en 
envergure  de  II  *  .84  a  T)  11  .99.  Nous  allons  considere  pour  un  meme  nombre  de  Mach  de  l'aile 
M  .  *  0.78  des  oscillations  instationnaires  autour  des  champs  stat ionnaires  qui  ont  d'ores  et 

dejSefe  calcuies  au  paragraphe  II.  Nous  pr^senterons  diverses  comparaisons  pour  bien  mettre  en 
Evidence  I'importance  de  la  correction  tr id imens ionne 1 le . 


III.l  -  Comparaisons  pour 


ai  le 


0.78 


ai  le 


aile 


0° 


Le  champ  de  pression  stationnaire  qui  sert  de  point  de  depart  &  Involution  instationnaire 
est  celui  prlsentl  en  figure  6.  La  frequence  d *  oscillation  de  la  gouverne  eat  de  20  Hz  (ce  qui 
correspond  A  une  frequence  rlduite  UC/l^c.24). 

L'amplitude  d 1 oscil lation  de  la  gouverne  est  1°.  Si  l'on  applique  directement  la 

mlthode  instationnaire  sans  introduire  les  corrections  tridimensionnelles  (nous  noterons  2.0  ) 

la  comparaison  des  premiers  harmoniques  calculus  et  mesurls  du  coefficient  de  pression  normalise 
A  l'extrados  (  Cp^  )  est  prlsentle  en  figure  8.  On  observe  que  le  module  calculi  est  netteraent 
surlvalul  (un  facteur  dix  dans  la  partie  amont  du  profil). 

Si  maintenant  nous  introduisons  les  corrections  tridimensionnelles  (nous  noterons  3D  ), 
pour  ce  faire  nous  utilisons  le  systlme  d'lquations  (34)  pour  dlfinir  le  pompage  et  le  tangage 
induits  qui  se  superposent  au  mouvement  de  la  gouverne,  nous  obtenons  la  comparaison  de  la  figure  9. 
Le  module  calculi  e9t  fortement  rlduit  par  cette  correction  mais  il  reste  cependant  trop  fort. 

Si  on  considlre  les  valeurs  obtenues  par  la  mlthode  linlaire  des  doublets  dont  le  rlsultat  est 
Igalement  portl  en  figure  9  on  arrive  A  la  meme  conclusion.  On  peut  done  penser  qu'on  est  1&  en 
presence  d'effets  visqueux  importants  qui  comme  on  le  sait  ont  tendance  A  rlduire  l'efficacitl 
instationnaire  sur  la  gouverne.  A  partir  des  travaux  de  GRENON  et  al.  [4]  sur  les  possibility 
de  couplage  du  fluide  parfait  et  de  la  couche  liraite  on  peut  essayer  de  prendre  en  compte  une 
partie  des  effets  visqueux  en  consid^rant  que  l'excitation  efficace  n’est  qu'une  certaine  fraction 
de  l'excitation  de  la  gouverne.  En  effet,  au  raoins  en  subsonique,  le  principal  effet  de  la  couche 
limite  instationnaire  est  de  rlduire  le  module  du  coefficient  de  pression  instationnaire  sans  trop 
affecter  la  phase.  En  figure  10  nous  prlsentons  le  rlsultat  d'un  calcul  (P.P.T  3D  +  "CL")  obtenu 
par  une  reduction  arbitraire  du  dlbattement  de  la  gouverne  de  25  %  et  sur  l'extrados  du  profil 
seulement.  Le  choix  d'un  abattement  sur  l'extrados  seul  vient  d'une  remarque  faite  h  partir  des 
essais.  En  effet,  dans  cette  configuration  le  module  mesurl  de  la  pression  instationnaire  est 
nettement  plus  faible  sur  l'extrados  que  sur  l'intrados  ce  qui  nous  a  serobll  £tre  le  fait  d'une 
predominance  des  effets  visqueux  du  cote  extrados .  L' analyse  du  resultat  obtenu  montre  bien  que  le 
module  de  la  pression  instationnaire  est  rdduit,  mais  de  manifere  insuf f isante .  Le  comportement 
de  la  phase  va  lien  dans  le  sens  de  ce  qui  est  observe  experimentalement  mais  la  correction  est  trop 
importante  principalement  au  bord  d'attaque.  On  peut  done  penser  qu'un  abattement  sur  les  deux 
faces  du  profil  et  different  de  part  et  d'autre  aurait  permis  d'obtenir  une  bonne  comparaison 
theorie-experience .  Il  semble  done  bien  qu'une  meilleure  prediction  des  champs  de  pressions 
instat ionnaires  passe  par  la  prise  en  compte  de  la  couche  limite  instationnaire. 

Si  l'on  considere  les  valeurs  de  la  portance  et  du  moment  quart  avant  pour  ce  cas  on 
obtient  le  tableau  suivant  : 


Ce* 

C-m6 

Module 

Phase 

Module 

Phase 

P.P.T.  2D 

3.85 

-  25° 

1.05 

+  3° 

P.P.T.  3D 

1.76 

-  4° 

.71 

♦  5° 

ESSAI 

1.04 

5° 

.45 

*  10 

DOUBLETS 

1.77 

-  1° 

.72 

4° 

P.P.T. 

3D  +  "C.L" 

1.27 

4° 

.57 

6° 

Plus  que  les  valeurs  locales  deB pressions  sur  le  profil  ce  sont  ces  grandeurs  qui  ont  de 
1' importance  si  l'on  s'intdresse  par  exemple  aux  calculs  de  flottement.  On  remarquera  tout  d ' abord 
que  si  l'on  ne  tient  pas  compte  des  effets  tridimensionnels  (P.P.T.  2D)  le  resultat  est  tr£s  loin 
des  essais .  La  prise  en  compte  des  effets  tridimensionnels  (P.P.T  3D)  amlliore  sensiblement  la 
comparaison  mais  ne  fournit  pas  une  meilleure  approximation  qu'un  programme  tridimensionnel  linlaire 
(ici  les  doublets).  Par  contre  la  prise  en  compte  d'effets  visqueux  (P.P.T  3D  ♦  "CL")  bien  que 
trls  grossi&rement  Ibauchle  ici  fournit  une  correction  appreciable  de  la  portance  et  du  moment.  On 
doit  done  accentuer  notre  effort  pour  un  couplage  du  calcul  de  fluide  parfait  et  de  la  couche 
limite  en  instationnaire.  Il  est  bien  sur  Evident  qu'un  tel  couplage  ne  peut  pas  se  faire  & 
partir  d'une  mlthode  linlaire  de  type  doublets,  ou  les  effets  d'lpaisseur  sont  nlgligls  et  oCj 
l’on  ne  tient  pas  comple  des  effets  non  linlaires  comme  les  ondes  de  choc  qui  interagissent  fortement 
avec  la  couche  limite . 

III. 2  -  Comparaisons  pour  Maile*  0.78,  0»aile  «  2°  ,6ail^  ~  6° 

Le  champ  stationnaire  correspondent  h  cette  configuration  est  celui  de  la  figure  7.  La 
frequence  d ' osci 1 lat ion  de  la  gouverne  est  ici  de  10  Hz.  L'amplitude  d ' osci llation  de  la  gouverne 
est  6 |  ■  1.9°.  Nous  prlsentons  en  figure  11  le  rlsultat  obtenu  par  le  calcul  dans  le  cas  ou 
l'on  ne  tient  pas  compte  des  effets  tridimensionnels  instationnaires .  Ici  encore  les  modules  du 
coefficient  de  pression  sont  trls  nettement  surlvaluls.  On  peut  noter  que  le  module  du  premier 
harmonique  du  coefficient  de  pression  est  trks  important  aux  environs  de  25  Z  de  la  corde  (e'est 
la  zone  qui  est  balayle  par  l'onde  de  choc),  le  module  dlborde  largement  du  cadre  du  dessin 
puisque  le  module  calculi  atteind  la  valeur  27. Cette  zone  est  Igalement  le  silge  des  principales 
non  linlaritls,  ainsi  localeraent  le  module  du  second  harmonique  du  coefficient  de  pression  peut 
valoir  jusqu'S  14,  et  le  troisilme  jusqu'd  8.  Sur  les  valeurs  intlgrles  de  portance  et  de  moment 
(au  quart  avant)  le  module  du  second  harmonique  reprlsente  cependant  moins  de  1  Z,  du  module  du 
premier  harmonique. 
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Si  l'on  applique  la  correction  tridimensionnel le  instat ionnai re  du  paragraphe  II. 3  on 
obtient  le  r^aultat  presente  en  figure  12.  L'excursion  du  choc  est  ici  sensiblement  r^duite  et  si 
l'on  observe  un  pic  de  pression  autour  de  25  X  de  la  corde  la  largeur  de  ce  pic  est  sensiblement 
plus  faible  que  dans  le  cas  precedent  figure  11.  On  remarquera  que  le  module  du  pic  de  pression 
calculi  est  en  bon  accord  avec  celui  observe  exper imentalement .  Paradoxalement ,  les  effete  visqueux 
semblent  ici  etre  sensiblement  moins  importance  puisque  les  modules  calculus  sont  sensiblement  en 
accord  avec  ceux  mesur^s.  On  doit  probablement  cela  au  bvaquage  moyen  de  la  gouverne  (  6  A{\e  *  ”  6°) 

qui  equipe  I'aile,  ce  braquage  est  favorable  k  la  couche  limite  de  l'extrados  au  moins  au  niveau  de 
la  gouverne.  En  ce  qui  concerne  la  phase  1‘ accord  est  relativement  bon  sau/  dans  la  region  du  choc, 
on  peut  voir  Ik  <?galement  un  effet  visqueux.  En  effet  la  phase  exper imentale  evolue  qual itat ivement 
comme  la  phase  calcuiee  mais  il  serable  qu'elle  ne  soit  pas  capable  d'evoluer  aussi  rapidement,  ce 
lissage  peut  etre  corapris  conine  le  r^sultat  de  1  *  interact  ion  de  I'onde  de  choc  et  de  la  couche 
limite  qui  n'est  bien  sur  pas  pris  en  compte  dans  le  calcul. 

IV  -  POSSI BILITES  OFFERTES  POUR  UN  CALCUL  DE  PREDICTION  DU  FLOTTEMENT 

La  prevention  d'un  eventuel  flottement  sur  les  ailes  supercritiques  qui  djuipent  les  avions 
de  transport  actuela  s'appuie  encore  aujourd'hui  sur  des  m£thodes  lin^aires  comme  la  mlthode  des 
doublets.  Ces  m^thodes  qui  ne  permettent  pas  de  tenir  compte  des  ondes  de  choc  sont  done  sujettes 
A  caution  dans  le  domaine  transsonique .  II  nous  semble  done  n£cessaire  d'aller  plus  loin  dans  la 
prediction  en  utilisant  des  m^thodes  non  lin^aires  du  type  de  celle  qui  a  ete  presentee  ici,  soit 
pour  justifier  1 'ut i 1 isat ion  des  methodes  lin^aires,  soit  pour  faire  ressortir  leurs  insuf f isances . 

Un  tel  calcul  s'il  apparait  de  prime  abord  comme  particul ifrrement  complexe  peut  beneficier 
d'un  certain  nombre  de  remarques  : 

-  linearite  et  harmonicite 

En  transsonique  comme  nous  venons  de  le  voir  tout  part icu 1 i^rement  au  passage  des 
ondes  de  choc  les  coefficients  de  pression  ont  une  Evolution  qui  n'est  ni  lineaire  ni  harmonique. 
Cependant  les  coefficients  globaux  de  portances  et  moments  peuvent  etre  considers  comme  harmoniques 
et  lin^aires. 

Cette  remarque  qui  s'appuie  sur  les  r^sultats  des  calculs  pr^sent^a  ici  ou  en  reference 
^5,6]  rejoint  des  observations  exp^r imentales  dont  cellos  de  la  re f<*renceQl~J 

-  superposition 

Le  principe  de  superposition  des  mouvoments  estegalement  fondamental  en  aeroeiast ic ite , 
il  permet  de  limiter  le  nombre  dea  calculs  en  ne  considerant  que  des  modes  e lementai res(un  pompage, 
un  tangage  et  une  gouverne).  Par  superposition  on  peut  &  partir  de  ces  modes  simples  reconstruire 
un  autre  mode.  Par  exemple  un  tangage  autour  de  la  mi-corde  superpost*  k  un  pompage  donnera  les 
memes  forces  ins  tat ionnaires  qu'un  tangage  autour  du  quart  avant .  Dans  le  cas  de  profils  non 
decolies  ce  principe  reste  valable  et  a  ete  verifie  exp^r imentalement  en  reference  £l2j  . 

On  peut  done  envisager  un  ( alcul  de  flottement  en  regime  transsonique  de  manure 
sensiblement  plus  simple.  A  partir  du  calcul  non  lineaire  dans  un  certain  nombre  de  tranches 
representatives  de  I'aile  on  determine  les  coefficients  de  portances  et  moments  lies  aux  modes 
eiementaires .  Puisque  ces  coefficients  varient  1 ineairement ,  que  seul  le  premier  harmonique  est 
significatif  et  que  le  principe  de  superposition  reste  applicable  on  peut  construire  la  matrice 
des  coefficients  aerodynamiques .  On  peut  alors  resoudre  le  syst^me  lineaire  traditionnel  qui 
d^coule  de  liquation  du  flottement  ou  interviennent  outre  la  matrice  des  coefficients  aerodynamiques 
les  matrices  de  masse,  de  raideur  et  eventuel lement  d ' amort i ssement  ;  A  ce  niveau  on  n'a  done 
plus  k  se  pr^occuper  des  phenomfcnes  non  lin^aires. 

Pour  1' instant  un  tel  calcul  est  en  cours  de  developpement  k  la  division  de  1 ' aeroeiast i- 
cite  de  l'ONERA.  Dans  un  premier  temps  ne  disposant  pas  encore  d'un  calcul  couple  avec  la  couche 
limite  instationnaire  la  matrice  des  coefficients  aerodynamiques  est  construite  sans  tenir  compte 
des  effets  visqueux.  Cette  hypothfcse  nous  semble  acceptable  bien  que  nous  ayons  vu  au  paragraphe 
Ilt.l  que  les  effets  visqueux  instat ionnaires  etaient  importants,  en  effet  dans  un  calcul  de 
flottement  (sur  des  ailes  de  ce  type)  la  participation  des  gouvernes  est  trfcs  faible  et  l'on  peut 
penser  que  les  modes  £l£mentaires  sont  moins  sensibles  aux  effets  visqueux  qu'un  mode  de  gouverne. 

V  -  CONCLUSION 

L'lvaluation  des  champs  de  pression  instationnaire  induit  sur  I'aile  ZKP  par  1 'oscilla¬ 
tion  de  la  gouverne  a  ete  entreprise  k  l'aide  du  programme  bidimensionnel  de  resolution  de  1 'equation 
non  lineaire  des  petites  perturbations  transsoniques .  Apr&s  avoir  introduit  certaines  corrections 
propres  A  generer  un  champ  de  pression  stationnaire  correct  la  notion  d'angle  induit  tridimensionnel 
instationnaire  a  ete  developpe.  Cet  angle  induit  evalue  A  partir  d'un  calcul  tridimensionnel  de 
doublets  a  ete  introduit  dans  la  methode  non  lineaire  pour  lui  apporter  une  correction  tridimension¬ 
nel  le.  Cette  correction  s'est  averAe  tr£s  intereasante  meme  dans  le  cas  ou  l'on  se  trouve  en 
presence  d'ondes  de  choc,  en  effet  si  l'on  ne  tient  pas  compte  des  effets  tr idimensionnels 
instationnaires  les  pressions  calcuiees  sont  trfcs  sensiblement  surAvaluees  par  rapport  aux  essais. 

L* importance  des  effets  visqueix  a  egalement  ete  soul igneelors  de  la  comparaison  aux  essais. 
11  semble  souhaitable  de  disposer  rapidement  d'un  calcul  couple  entre  la  couche  limite  instationnaire 
et  le  calcul  non  lineaire  pour  obtenir  une  meilleure  prediction. 


Les  possibility s  offertes  par  cette  approche  en  regard  de  la  prediction  du  flottement  sont 
ygalement  analyses.  11  semble  qu'un  tel  calcul  soit  envisageable  puisque  les  portances  et  moments 
satisfont  au  principe  de  superposition  et  gardent  un  compor tement  harmonique  et  lin£aire.  Ce  calcul 
actuellement  en  cours  &  la  Division  de  1  * Aero£ las t ic i t£  de  I'ONERA  permettra  d'apporter  une  r^pcnse 
au  probl$me  que  pose  actuellement  1 'utilisation  des  m^thodes  linlaires  sur  les  ailes  supercritiques 
Iquipant  les  avions  de  transport  actuels. 
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Fig.  1  —  Evolution  du  vriltage  gEomEtrtque  et  du  vnllage 
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me  methoue  ues  tranches  qui  prend  en  Applications  de  la  proposent  une  mdthode  des  tranches  qui  prend  en  Applications  de  I; 
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